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Abstract 

In this paper, we provide two approximations in law of operator fractional Brownian motions. 
One is constructed by Poisson processes, and the other generalizes a result of Taqqu (1975). 
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1. Introduction 

Self-similar processes, first studied rigorously by Lamperti [TB] under the name "semi-stable" , 
are stochastic processes that are invariant in distribution under suitable scaling of time and space. 
There has been an extensive literature on self-similar processes. We refer to Vervaat [3D] for general 
properties, to Samorodnitsky and Taqqu [28] [Chaps. 7 and 8] for studies on Gaussian and stable 
self-similar processes and random fields. 

The definition of self-similarity has been extended to allow for scaling by linear operators on K d , 
and the corresponding processes are called operator self-similar (o.s.s) processes in the literature. 
See Laha and Rohatgi [T!5] , Hudson and Mason [T7] , and Sato [5DJ . Various examples of operator 
self-similar Gaussian and non-Gaussian processes have been constructed and studied by Maejima 
and Mason [21] and Mason and Xiao [22]. The theory of operator self-similarity runs somewhat 
parallel to that of operator stable measures and is also related to that of operator scaling random 
fields. See, for instance, Meerschaert and Schefiler [51], Bierme et al.[2] and the references therein. 

Let End(M. d ) be the set of linear operators on M. d (endomorphisms) and let Aut(M. d ) be the set of 
invertible linear operators (automorphisms) in End(M. d ). For convenience, we will not distinguish 
an operator D 6 End(M. d ) from its associated matrix relative to the standard basis of IR d . 

In this paper, we will use the following definition of o.s.s. processes, which corresponds to that 
of Sato [26] , but is stronger than that of Hudson and Mason [17] 

An R d — valued stochastic process Y = {Y(t}} is said to be operator self-similar if it is stochas- 
tically continuous, and there exists a D G End{R d } such that for every c > 

Y(ct) = c D Y(t), t <= R, (1.1) 

where = denotes equality of all finite-dimensional distributions, and 

c D = exp((\ogc)D)) = ^(logc) k D k . 

fc=0 
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Any matrix for which (11.11) holds is called an exponent of the o.s.s process Y. 

We say that a process Y = {Y(t)} has stationary increments (s.i.) if for every b € K 

Y(-+b)-Y(b)±Y(.)-Y(0). (1.2) 

One of examples of operator self-similar processes is the operator fractional Brownian motion 
(OFBM). OFBMs are mean-zero, o.s.s., Gaussian processes with stationary increments. They are 
of interest in several areas and for reasons similar to those in the univariate case. For example, they 
are obtained and used in the context of multivariate time series and long range dependence (see, for 
example, Chung [3] , Davidson and de Jong [6] , Davidson and Hashimzade [7] , Dolado and Marmol 
|12) . Robinson [25], and Marinucci and Robinson 23 J. They are also studied in problems related 
to, for example, queuing systems and large deviations (see Delgado jS], and Konstantopoulos and 
Lin[20 ). In particular, Mason and Xiao [52] studied the sample properties of a particular class 
of OFBMs. Didier and Pipiras [TO] [11] studied the basic properties of OFBMs, such as the time 
reversibility, the behavior of the spectral density around zero, and so on. 

On the other hand, weak convergence to FBM processes has been studied extensively since the 
works of Davydov 5 and Taqqu [29) . In recent years many new results on approximations of 
FBM processes have been established. For example, Enriquez [13] showed that a FBM can be 
approximated in law by appropriately normalized correlated random walks. Delgado and Jolis [8] 
proved that the law of a FBM can be weakly approximated by the law of some processes constructed 
from a standard Poisson process. 

Despite a growing interest in OFBMs, there is little work that studies weak limit theorems 
for OFBMs. In this paper we will study weak limit theorems for OFBMs. We will extend two 
approximation results for the one-dimensional FBM to the multivariate case of OFBM. The rest 
of this paper is organized as follows. In Section 2, we recall some preliminaries and present the 
main results of this paper. Based on Poisson processes, we prove weak limit theorem for OFBMs 
in Section 3. Based on a stationary sequence, we study weak convergence to OFBMs in Section 4. 



2. Preliminaries 

We first recall some facts we need later. Throughout this paper, we will use \\x\\e to denote 
the usual Euclidean norm of X € R d . For A e End(R d ), let \\A\\ = maxy !1 .|| JS= i denote the 

operator norm of A. It is easy to see that for A,Be End(R d ), 

||AB|| < ||A|| • ||S|| , (2.1) 

and for every A = (Ay) dxd € End(R d ), 

max JAij] < \\A\\ < rfi max |A y |. (2.2) 

Let cr(A) be the collection of all eigenvalues of A. We denote 

X A = min{i?eA : A € cr(A)} and = max{i?eA : A E er(A)}. (2.3) 

As standard for the multivariate context, we furthermore assume that the OFBM is proper. A 
random variable in M. d is said to be proper if the support of its distribution is not contained in a 
proper hyperplane of M. d . 

Let x' denote the transpose vector of a; € K d , and B* be the adjoint operator of B. Let 
x + = max{x, 0} and X- — max{— x, 0}. 

Didier and Pipiras 10 established the integral representations of OFBMs in the spectral domain. 
The following lemma comes from Didier and Pipiras [10] [Theorem 3.1]. 
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Lemma 2.1 Let D be a linear operator on ~R d with < Ad,A,d < 1. Let X — {X(t)} be an 
OFBM with o.s.s. exponent D. Then X admits the integral representation 

X(t) = f flzl ( x-P-Va + xZ (D ~ i} A] W{dx) (2.4) 
J M ix V ) 

for some linear operator A on C d . Here, A denotes the complex conjugate and 

W{x) := W 1 (x)+iW 2 (x) 

denotes a complex-valued multivariate Brownian motion such that W\{— x) = Wi(x) and W 2 {—x) = 
— W 2 (x), W\{x) and W%{x) are independent, and the induced random measure W(x) satisfies 



W(dx)W*(dx) 



= dx, 



where W* is the adjoint operator of W . 
Remark 2.1 For fixed t € K, let 



F(x)= e ^l(x- + ^A + xZ {D - i) A 



ix 

= F 1 (x)+ iF 2 (x). 

Since F(x) = F(-x), F^x) = Fi(-x) and F 2 (x) = -F 2 (-x). 

Inspired by Samorodnitsky and Taqqu [28] [Chap. 7], up to a multiplicative constant, we can 
rewrite {X(t)} as follows. 

/•OO /"OO 

X(t)= I G 1 (x,t)W 1 (dx)+ I G 2 {x,t)W 2 (dx), (2.5) 

where 



o Jo 



and 



„ SmtX _(r>_L) . costx - L _m_I> . 

X X 

sinte _(r>-L\ . 1-costx . 

X X 

A = A 1 +iA 2 . 



Remark 2.2 Using the same method as Mason and Xiao [55], we can get that for i = 1,2, 

/ ||Gj(u, t) \\ 2 du < 00. 

In this paper, we study weak limit theorems for the OFBM X. We hrst recall some results about 
weak convergence to FBMs. Stroock [27] showed that 

Lemma 2.2 Let {N n (t), n = 1, 2 • • • } be a sequence of Poisson processes with intensity n, and 
construct the continuous processes {Y„ (t)} by 

Y n (t) =y/n f (-l) N ^ s) ds, < t < 1. (2.6) 

Then the laws of {Y n } converge weakly, in C([0, 1]), to the law of a Brownian motion, as n — > 00. 
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Using the same method as Delgado and Jolis [5], one can easily get 

Corollary 2.1 Let N n (t) and N n (t) be two independent Poisson processes with intensity n, and 
construct the processes {Y n (t)} by 

Y n (t) = Vn-[ ^^(-l) N ^d X + ^[ (-!)*.<•>*,, (2-7) 

JR + x ' JR + x 1 

where H G (0, 1). Then the laws of {Y n (t), n — 1, 2, • • • }t e ro, i] converge weakly in C([0, 1]), to the 
law of a fractional Brownian motion with index H G (0, 1), as n — > oo. 

Inspired by Corollary 12.11 we want to show that the OFBM X given by (|2.5p can also be 
approximated by a sequence of processes similar to (|2.7|) . 
Let 9° n {t) = vM-l)^* 5 and 

=(#(*),... >e d (t))', (2-8) 

where 6 l n {t), i G {1, • • • , d} are independent copies of Furthermore we define 

] n ( S )d s =(J*9 1 n ( s )ds,--- ,Je d n {s)ds)'. (2.9) 

Inspired by Corollary 12.11 and the construction of {Y n }, we define the sequence {X n (t)} n( =m as 

X n (t)= G 1 (x,t)9 n (x)dx + G 2 {x,t)9 n (x)dx, (2.10) 
Jo Jo 

where 9 n (x) is an independent copy of 9 n (x). Then 

Theorem 2.1 The laws of {X n (t), n — 1, 2, • • • }t£[o, i] given by 12. 1 0\) in C d [0, 1] converge weakly 
to the law of the OFBM X given by pO]) , osn->oo, where C d [0, 1] = C([0, 1], R d ) . 

On the other hand, Taqqu |29j showed that a FBM can be approximated in law by normalized 
partial sums of stationary random variables. Since OFBMs are multivariate extensions of FBMs, 
it is interesting to extend Taqqu's [29] result to the OFBM case. 

Before we state our result, we first introduce the following notation. 



and 



{A(n)} = {{A t3 (n)) dxd } e End(R d ) 
{B(n)} = {(%(n)) dxd } G End(R d ) 



are asymptotically equivalent, as n —¥ oo, if for any i,j G {1, • • • one of the following cases 
holds: 

(i) There exists N e N such that for all n> N , 

Bij{n) ^ and lim Aj 7 (n)/Sj,(n) = 1. 

n— too 

(ii) There exists N% G N such that for all n > N% , 

Bij(n) = and Aij(n) = 0. 
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We denote this as A(n) ~ B(n), as n — > oo. 

Recall that a process Y = {Y(i)} is time reversible if it satisfies 

Y{t) = Y(-t), t G R. (2.11) 

The author refers the readers to Didier and Pipiras [TU] for time reversible Gaussian processes with 
stationary increments. 

Due to the following lemma, we only can extend Taqqu's result to the time reversible case. 

Lemma 2.3 Let {Zi,i — 1,2,---} be a stationary proper mean-zero Gaussian sequence ofM. d - 
valued vectors, and 

\nt\ 

1=1 

If all the finite- dimensional distributions of{S n (t)} converge to the corresponding finite- dimensional 
distributions of a mean-zero proper Gaussian process Z = {Z(t)}. Then Z is time reversible. 

Proof: Define 

E[Z(t)Z'(s)]=R z (t,s). (2.12) 

Since 

S n (t) Z(t),as n — ^ oo, 

F D 

where denotes convergence of the finite-dimensional distributions, we get that 

lim E[S n (t)S' n (s)] =R z (t,s). (2.13) 

71— >OC 

Since {Zi} is a stationary Gaussian sequence, we can get 

E[S n {t)S' n (s)] =E[S n (s)S' n (t)] (2.14) 
It follows from (|2~T2]) and (|2~T4|) that 

R z (t,s) = R z (s,t). (2.15) 

Since {Zi} is a stationary Gaussian sequence, and Z is a Gaussian process, we get that Z has 
stationary increments. Therefore, by Proposition 5.1 in Didier and Pipiras [10] and (|2.15[) . we get 
that the lemma holds. □ 

Remark 2.3 It follows from Proposition 5.1 in Didier and Pipiras [10] that if the OFBM X is 
time reversible, then 

E[x(t)x'( s )] = \ [\t\ D r\t\ D ' + \ s \ D r\s\ D * - \t - s\ D r\t - S \ D ' 

Now we state our result as follows. 

Theorem 2.2 Let {Zi,i = 1,2, •••} be a stationary proper mean-zero Gaussian sequence ofR d - 
valued vectors. We define 

r(i,j)=E[Z i Z^] = (r kq (\i-j\j) 



ixd 



Suppose that 

N N 



/Z/Z^'ri ~ KBN D TN D * B*, asN^oo, (2.16) 
i=i i=i 
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where T = ~E[X(1)X '(!)], B G Aut(R d ) and K > is a positive number. Then 



Q N {t) = d N 



[Nil 

E 



(2.17) 



with djy ~ CN B , converges weakly, as N — > oo in T) d [0, 1], up to a multiplicative matrix from 
the left, to the time reversible OFBM X given by with A 2 A\ = A X A* 2 , where C G Aut(R d ), 

and V d [0,l] = V([0, l],TR d ). 



We will prove Theorem 12.21 in Section 4 . 

To end this section, we give a technical lemma which comes from Maejima and Mason [3T] . 

Lemma 2.4 Let D G End(R d ). If > and r > 0, then for any 5 > 0, there exist positive 
constants K\ and K 2 such that 



I dm 1 Kir x °~ s , for allr < 1, 
H - 1 K 2 r AD+s , for allr > 1. 



(2.18) 



In the rest of this paper, most of the estimates contain unspecified constants. An unspecified 
positive and hnite constant will be denoted by K, which may not be the same in each occurrence. 
Sometimes we shall emphasize the dependence of these constants upon parameters. 

3. Weak limit theorem based on Poisson processes 



In this section, we prove Theorem 12.11 We first state some technical results. In order to prove 
the main result, we need a tightness criterion for vector-valued stochastic processes. 

Lemma 3.1 Let |Z„(i)} n( _ N be a sequence of stochastic processes inT> d [0, 1] satisfying: 
(i) for every n G N, Z n (0) = a.s.; 

(ii) there exist constants K > 0, /3 > 0, a > I and an integer Nq G N such that 

0' 



E 



Z n (t) - Z n (s) 



< K(t-s) a ,n> N a . 



(3.1) 



Then {Z n (t)} is tight in V d [0, 1]. 



The proof of Lemma 13.11 is classical (see Billingsley |T| [Chap. 3] and Ethier and Kurtz [Tl] 
[Chap.3]). Here we omit the proof. 

Next, we show that {X n (t)} is tight in C d [0, 1]. We need the following lemmas. 

Lemma 3.2 For any even m G N, there exists a constant K(m,d) > such that for any measur- 
able function f : M + -> End(R d ) with J R \\f(u)\\ 2 du < oo, 



E 



f{u)6 n (u)du <K(m,d) 



\f(u)\\ 2 du 



(3.2) 



where n (-) is given by 

Proof: Let f(u) = (&(«)) rfxd = (/H«), ■ ■ • . f d (u))\ where P{u) = (f a (u), ■■■ , f ld {u)). Then 



f(u)6 n (u)du = ( [ f(u)e n (u)du, ■ ■ ■ , f f d (u)6 n (u)du)' 

V ./"IQ> , ./tip . / 



(3.3) 
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It follows from (13.31) that 



a - 

<K(m)J2^[\ / f{u)6 n {u)d 



E 



E / f{u)0 n (u)du 
For each i G {1, • • • , d}, we have 

f(u)9 n (u)du m ] <X(to)Ve[| / f i:j (u)ei(u)di 

Since J R | |/(-u)| | 2 dw < oo, it follows from ([272]) that 

|^-(u)| 2 du< / ||/( U )|| 2 d U <co. 
Using the same method as in Delgado and Jolis [8], we have 

/y(u)0£(u)du]l m < Jf(m)( / |/«(u)| 2 du 



(3.4) 



(3.5) 



(3.6) 



(3.7) 



□ 

Lemma 3.3 Choose 5 > 0, suc/i i/iai A^j — 5 > and +5 < 1, then for any < s < t < 1 and 
even m £ N, we have 



Combining (j2~2j) . (j3~4l) . (f3~5|) and ([3~7| . we get that the lemma holds 



E 



X n (t) - X n (s) 



< K(m,5, d)(t-s) 



mH 



(3.8) 



where H = Xjj — 5. 



Proof: Using the inequality (a + b) m < 2 m 1 (a m + b m ) for a, b > 0, we get that there exists a 
constant K(m) > such that 



E 



X n (t) - X n (s) 



< K{m)l 



(G 1 (x,t)-G 1 (x,s)y n (x)da 



(G2 (x, t) — G 2 (x,s)) n (x)dx 



where 



and 



K{m)Ii{m) + K(m)l2(m), 



(G 1 (x,t)-G 1 (x,s)y n (x)da 



(G 2 (x,t) - Gz(x, s)\d n {x)dx 



(3.9) 



h(m) = E 



h{m) =E 



Since 9 n {x) is an independent copy of 9 n (x), we only deal with Ii(m). I2(m) can be done in the 
same method. 

By Lemma 13.21 and Remark 12.21 we get that 



h 



(to) < K(m,d)( j \\Gi{x,t) - Gi(a;,s)|| 2 dx) 



(3.10) 



On the other hand, 



\G\{x,t) — Gi(x, s)|| 2 = (sinta — sin sx) (-) D+ ^ 

+ (costa — cos sx)( — ) D+ i 

1 i 2 
<K (sinte — sinsx)(-) D+5 

v ' x 

+K (costx — cossx)(-) D+ ^ 



By (13TU1) and (IXTT1) . 

Ii(m) < K(m) 
+K(m) 

For the sake of conciseness, we define 



I / . , . w 1 x D+i II 2 

(sin te — sin sx)( — ) 2 



| (cos tx — cos sa; ) ( — ) 
x 



dx 



dx 



and 



U 2 {m) = 

It follows from Lemma \2. 41 that 



(sinte — sinsx)(— ) 
x 



(coste — cossx)( — ) 
1 x 



dx 



dx 



U 2 (m) <K(m,S,d)( I (costx~cossx) 2 (-) 2(AD+5+i) 



d/x 



+K(m,6,d)(j (costx-cossx) 2 (-) 2{XD S+ ^dxJ 
The first term on the r.h.s. of (|3.13[) can be bounded from above by 



K(m,5,d)( f sm 2 ((t - s)x/2)(-) 2{AD+S+ ^dx 
* J o x 



< K(m,S, d)(t - s) m< - AD+S) 

< K{m,S, d)(t ~ s) m(AD+s *>( 



00 sin 2 x 



2(A D +5+i) 



dx 



o x 



1 • 2 

sin x 



o a; 



2(A£,+<5+i) 



rzdx + 



1 



1 x 



2(A D +<5+±) 



d.r 



The second term on the r.h.s. of (|3.13|) can be bounded from above by 



1n2(A d -5+!) 



K(m,S,d)( / sin 2 ((<-s)a;/2)(-) 

< K(m,5, d)(t-s) m ^ D - s ^ 

< K(m,S, d)(t- sy^n-s)^ 



dx 



■ 2 

sin x 



2(A D -<5+i) 



dx 



1 • 2 

sin x 



o x 



rrdx + 



1 



1 X 



2(A D -5+i) 
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It follows from (|3~T5j) to (157151) that 

U 2 (m) < K(m,6,d){t- s) mH , 

where H = \d — 5. 

Using the same method as the proof of (j3.16j) , we get that 

Ux(m) < K(m,6,d)(t- s) mH . 
It follows from (|3T2j) . (|3~T6l) and (13T7)) that 

him) < K(m,6,d)(t- s) mIi . 
Using the same method as the proof of (J37T8J) , we get 

J a (m) < K(m,S,d)(t- s) mH . 
The lemma follows from ([31?]) , (fXT5|) and (gTTJJ) . 



(3.16) 

(3.17) 
(3.18) 

(3.19) 

□ 



Remark 3.1 It follows from Ethier and Kurtz [T3] [Chap. 3 Proposition 10.3] that Lemma I3~3l also 
shows that P(X n e C d [0, 1]) = 1. 

In the following, we will prove that {X n (t), t g [0, 1]} given by (12.10)) converges weakly in 
C d [0, 1] to the OFBM X given by (]2.5p . as n — > oo. In order to obtain it, we need the following 
lemma. 

Lemma 3.4 The laws of { J* 9 n {u)du} converge weakly in C d [0, 1] to the law of an M. d -valued 
Brownian motion, as n — > oo. 



Proof: For every even m € N and any < s < t < 1, 



E 



i(u)du 



i(u)du 



By Lemma 2 in Delgado and Jolis [S], we get 



E 



i(u)du ■ 



i{u)du 



< K(m)(t - s) ! 



(3.20) 



(3.21) 



By Lemma |3~T1 and Proposition 10.3 in Ethier and Kurtz P~4] [Chap. 3], (13.21)) implies the tightness 
01 {fo°n(u)du}. 

On the other hand, the convergence of all hnitc-dimcnsional distributions follows from Lemma [2.2 
and the well-known Cramer- Wold device (Jurek and Mason [16] [Chap.l] and Whitt [31] [Chap. 3]). 

Finally tightness plus convergence of the finite-dimensional distributions implies the weak con- 
vergence. □ 

The proof of Theorem 12. II is as follows. 

Proof of Theorem \2.1\ We first prove the tightness of the laws of {X n (t), t e [0, 1]}. From 
Lemma T3.31 for every even m € N, which satisfies mH > 1, inequality (|3.8j) holds. It follows from 
Lemma T3. II that {X n } is tight. 

Now we proceed with the identification of the limit law. It is sufficient to prove that for any 
g6N,Oi,-",o,6l and t\ ■ ■ ■ ,t q G [0, 1], 



(3.22) 
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We define for any me {1, • • • , q}, 

G 1 {t m ,x)0 n (x)dx= (x^(t m ),--- ,X*(t m ) 
G 2 {t m ,x)9 n (x)dx = (xl{t m ), ■ ■ • ,**(t m ))' 



1 (t m ,x)dW 1 (x) = (x\t m ),--- ,X d (t m ))\ 



and 

' G 



+ 

G 2 (t m ,x)dW 2 (x)= (l 1 ^),--- ,X d (< m ))' 



By the Cramer- Wold device, in order to prove (I3.22p we only need to show that for any 
(b 1 , ■■■ , b d )' € M d 

d q d q 

J2 Y a™VXi(t m ) J2 Y anV&dm), ( 3 - 23 ) 

j=l m— 1 j — 1 m— 1 

and 

d q d q 

Y Y a m b>Xi(t m ) ^YY a m b>P{t m ), (3.24) 

j — l m—l j — 1 m— 1 

since 6*„(u) and 6* n (u) are independent. Now we only prove that (13.23|) holds. p.24|) can be done 
in the same way. 

In order to simplify the notation, let 

Wi(x) = {Wl{x) 7 Wl{x) 7 --- ,Wf(x))', 

and 

Till T* 1 — I - -.ill HP \ — / T • ■ { 1 I <Y 1 , 

Uxd 

where f J (u,x) = (/ji (u, x), ■ ■ ■ ,f j d (u,x)). 
Therefore, we can rewrite (|3.23|) as follows. 



f(u,x) = Gi(u,x) = (fji(u,x)) d . = (^{u.x), ■ ■ ■ ,f d (u,x)) , 



where 

d q 



H n = YY Va ™ l fji(tm,x)eUx)dx = f F t { 
1=1 ro =l JR+ Jk+ 



x)6 l n (x)dx, 



d q „ . 

H l = YY bJam fji(t m ,x)dWl{x)= Ft 

1=1 m=l Jr + • 7r + 



(*W(aO, 



and 

d 9 

Fj(x) = X! X! bJa mf 3 i{t m ,x). 

j — l m — l 
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In order to prove p.25[) . it is sufficient to prove that for any (el, 



E 



exp 



a 



k=l 



a 
fc=l 



as n — > oo. 

Since 9*, k = 1, • • • , d are mutually independent, and so are W*, k = !,■■■ ,d, 



(3.26) 



E 



exp 



fe=i 



E 



< 



d 
fe=l 

E 



fe=l 

d 

JJ E[exp [iflT*]] - [] E[exp [i£ff fe ] 
* =i fc=i 

[exp -E^exp^ff 1 

d d 

JJ E[exp - H E[exp ^ff fc ] 



fc=2 



k=2 



By induction on fc = 2, • • • , d, we get that in order to prove (I3.26[) it is sufficient to prove that for 
k £ {!,-•• ,d}, 



as n — >• oo. 

By (|2.2j) . we get that for every G {1, • • • , d}, 

d q 



(3.27) 



f F£(u)du < I ffk(u,t m )du < oo. 



(3.28) 



Therefore, the proof of p. 271) follows the lines of the proof of Theorem 2.1 in Dai and Li 0]. 
Combining (|3~2%t and (pT27]> . we get that (13331) holds. Similarly, we get that (|3~1H1) holds. 
By Theorem 7.8 in Ethier and Kurtz [14] [Chap. 3] and Remark 13.11 we get that Theorem 12.11 
holds. The proof is done. □ 

4. Weak convergence based on a stationary sequence 

In this section, we will prove Theorem 12.21 We first show that {Qjv(i) : t £ [0, 1]} is tight in 
T> d [0, 1]. Before we prove tightness, we give a technical lemma. 

Lemma 4.1 Let {A(n)},{B(n)},{C(n)} and {D(n)} £ End{R d }. If 

A(n) ~ C(n), as n — > oo, 

and 

B{n) ~ D(n), as n — > oo, 

then 

A(n)B(n) ~ C(ri)D(n), 

as n — > oo. 

One can easily get that the lemma holds. Here we omit the proof. 
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Lemma 4.2 For any < s < t < 1 and even m € N, </iere exists N x G N smc/i i/iai 



E 



Qjv(t) - Qjv(s) 



< if(t- s )w(Az ) -«) ! N > Ni, 



where S > mf/i Ac — <5 > 0, and Qpf{t) is given by 
Proof: For any s < t <G [0, 1], we have 

LJV(t- s )j 



E 



Q N (t)-Q N (s) 



m 


= E 




E 















m 




= E 




Q N {t - s) 













since {i^} is a stationary sequence. 
In order to simplify the notation, let 



VN(t-s)\ 

d,N ^ Zi = (Q [jv(t-s)J i ' ' ' ' Q[N(t-s)}Y- 

i=l 

Since (a + b) m < 2 m - 1 (a m + b m ) for a, b > 0, we get that 



(4.1) 



(4.2) 



E 



} N (t) - Qn(s) 



< K{m) E 



i=i 



Q 



lN(t-s)l 



Since {Zi} is an Revalued Gaussian sequence, Q 



lN(t-s)\ 



is Gaussian. So, we get that E 



(4.3) 



Q 



IN(t-s)] 



is proportional to E 



Q 



lN(t-s)} 



It follows from (l4~3t that 



E 



Q N {t) - Q N {s) 



m 


< k(e 




E 







} N {t) - Qn{s) 



(4.4) 



Note that for an Revalued random variable Q = (Q 1 ,--- ,Q d )\ E[||Q|||;] equals the sum of 
diagonal entries of the correlation matrix. So it follows from (|2.2p that 



E 



] N (t)-Q N {s) 



m 


< K 


E 


E 







[Q N (t) - Qn(s)] [Q N (t) - Qn(s)]' 



By (f2~T6j) and Proposition 2.2.2 in Meerschaert and Schemer [24], 

lim E[Q N (t) - Q N (s)][Q N {t) - Qn(s)]' = KC{t - s) D T(t - s) D " C* . 

N— >oo 

It follows from (|431> and ([LB]) that there exists N x e N such that for all N > N 1; 



E 



Q N (t)-Q N (s) 



m 






< K 


E 





(t - s ) D r(t - s ) 



D* 



since ||-|| is continuous. 

By p~7) and Lemma [23 we get that for all N > N x , 



E 



Q N (t)-Q N {s) 



D\\^2 



<K\\{t-s) 



(t-s) 



D' 



(4.5) 



(4.6) 



(4.7) 



(4.8) 
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The proof is completed. □ 
Now we prove Th.eorem l2.2l 

Proof of Theorem W3i It follows from Lemmas l3~T1 and l4~2l that {Qiv(t)} is tight. 
Now we prove the convergence of all finite-dimensional distributions. For any p € N and 
ti, ■ ■ ■ ,t p £ [0, 1], <3jv(ti), • • • , QN(tp) are jointly Gaussian, since {Zi} is a Gaussian sequence. 
Define: 

RN(U,tj) = E[Q N {U)Q' N {tj)]. 
In order to simplify the notation, we define 

s N (t) = J2 Z >- 

i=l 



J' 



Since {Zi} is a stationary Gaussian sequence, we have that if t j > t 

E[QN(Jt i )Q , N(fi)] = ^d N \E[S N (t i )S , N (t i )] 

+E[S N (t j )S' N (t j )\ - E[S N (t t - t 3 )S' N {U - t^] 



"JN • 



(4.9) 



By ^M, 



[N tj \ [N tj \ 

E[S N (t 3 )S' N (t 3 )} = r (*' fc ) - KBN D tfTtf*N D 'B* 



i=l k=l 



as N — > oo. 
So, we get 



(4.10) 



lim R N (ti,tj) = R(ti,U), 

iV— >oo 



(4.11) 



where R(U,tj) = f C 

It follows from Remark [2J5I that if X is a time reversible OFBM, then 



\u\ u t\u\ u -\u-tj\vr\u-tj\ 



c* 



E 



CX(t)(CX(s))' 



1 



-c 



\t\ D T\t\ D * + \s\ D T\s\ D ' - \t - s\ D T\t - a\ D ' 



C* 



(4.12) 



Since Qjv(ii), ■ • • , Qjv(ip) are jointly Gaussian, it follows from (14.111) and (|4.12p that 

(Qjv(ti),--- ,Q^(*p)) ^^(^(ti),.-. ,ci(g). 

By Theorem 7.8 in Ethier and Kurtz [14] [Chap. 3] , we get that the convergence of finite-dimensional 
distributions and tightness ensure weak convergence of {Qat(£)}- The proof is completed. □ 

At the end of this paper, we consider the OFBM given by Mason and Xiao [35]. Mason and 
Xiao [22] introduced the OFBM Y = {Y(t)} as follows. 

Definition 4.1 Let D be a linear operator on K d with < Xfj, A^> < 1. For any t € K, define 



Y(t) 



1 x D+I/2 „~ 



(1 - cos(tx)) {-) ^ 1 dW{x) + / sinte • (-) ^ 1/z dW(x) 



l%D+//2 



(4.13) 



where W is a vector- valued Gaussian measure, and W is an independent copy of W . 
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By m| and P~T5|) . we see that, up to a multiplicative constant, the OFBM Y and the OFBM 
X with A — I have equal distributions. 
Define 

Y n (t)= f (l-cos(tx))(-) D+I/2 d n (x)dx+ [ sintx- (-) D+I/2 9 n (x)dx, 
Jr + x Jr + x 

where 9 n (x) is given by (|2.8|l and n (x) is an independent copy of 9 n (x). We immediately have 

Corollary 4.1 77ie Zaiws of {Y n (t),t € [0, 1]} converge weakly to the law ofY given by {J^.lSfy in 
C d [0, 1], as n — y oo. 

By (f47T3|) . we get that 

E[Y(i)Y'( S )] - i [|t| D f \t\ D * + \s\ D f\s\ D ' -\t- s\ D f\t - s\ D '] , (4.14) 

where f = E[Y(1)Y'(1)] . One can get that E[Y(t)Y'(s)] = E[Y(s)Y'(t)] . Therefore, by Propo- 
sition 5.1 in Didier and Pipiras 10 , the OFBM Y given by (|4.13p is time reversible. Therefore, 
from Theorem 12.21 we get that 

Corollary 4.2 Let {Zi, i = 1, 2, • • • } be a stationary proper mean-zero Gaussian sequence of K d - 
valued vectors. We define 

r{i,j)=E[Z i Z' j \ = (r kq {\i-3\)) t 

Suppose that 



dxd 



N N 

ft ~ KBN D TN D *B\ asN^oo, (4.15) 

i=l j=l 

where B £ Aut(R d ) and K > is a positive number. Then 

[Nil 

Q N (t) = d N Y, %U (4-16) 

i=i 

with dpf ~ CN~ D B~ x , converges weakly as N — > oo in T> d [0, 1], up to a multiplicative matrix 
from the left, to the OFBM Y defined by \£Wj. 

In the following, we give a special case where conditions in Theorem 12.21 and Corollary 14.21 are 
satisfied. 

Corollary 4.3 Suppose that D e Aut(R d ) is a real diagonal matrix with positive diagonal entries 
Ai, • • ■ , Xd- Let {Zi, i = 1, 2, • • • } be a stationary proper mean-zero Gaussian sequence ofM. d -valued 
vectors. We define r{\i — j\) = E[ZjZ'j\. Suppose that 



(i) Case \ < \ D < 1: 

(ii) Case < A D < \: 
and 



r(j) - Kfj 2D - 21 as j -> oo with K > 0. (4.17) 
r(j) - KTj 2D - 21 as j -> oo with K < 0, (4.18) 



-(0) + 2^r(j)=O, (4.19) 

3=0 



where O denotes the zero matrix. 
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Then 



(t) = d N J2 Z h (4-20) 



t=l 



with dft ~ N D C 1 , converges weakly as N oo, up to a multiplicative constant, to the OFBM 
Y given by ft4.13\ ), where C G Aut(R d ) is a diagonal matrix with diagonal entries 



V|2Ai-l|2Ai y/\2\ d -l\2\ d 
Proof: We will proceed to prove this corollary in two steps. 

Step 1. We first deal with the case that | < Ad < 1. Since D is a diagonal matrix, we get that 
j2D-2i - g s ^jj a ,-ii a g OIia i matrix. By some simple calculations, we get that the diagonal entries of 
y u ZJ are 

■2Ai — 2 •'2\ r i — 2 

J , ■ ■ ■ ,3 

For every i g {1, ■ ■ • ,d}, we get from Feller [15] [p. 281] and Taqqu 29 that 

^fc 2 ^- 2 ~(2A i -l)- 1 i 2Ai - 1 , as j^oo, (4.21) 



fe=i 



since | < Au < Ad < 1. 
So, by (H2J), 



E^ 27 ~Coi 2 ^, as i^oo, (4.22) 
fc=i 

where Co € Aiti(]R d ) is a diagonal matrix with diagonal entries 

(2AX-1)- 1 ,--- ,(2A d -l)" 1 . 

On the other hand, we have 

N N JV-1 j 

i — 1 J — 1 J=l 2—1 

Since D is a diagonal matrix, T is a diagonal matrix. Hence, we get from (14.17[) . (|4.22l) and 
P~2"5|) that, as N -> oo, 

JV 

.Bat - 2ifC f ^ j 213 " 1 - 2XC CifiV 2 - D , (4.24) 

where Ci G Aut(K d ) is a diagonal matrix with diagonal entries 

(2A1)" 1 ,--- ) (2A d )" 1 . 

Since Co and C\ are all diagonal matrices, we choose a diagonal matrix C G ^4ui(R d ) such that 
C 2 = CqC\. Since C,N D and T are all diagonal matrices, 

£tv - 2,fiXW I) f N D *C*. (4.25) 

Using the same method as the proof of Theorem [221 we get that {Q./v(i)} converges weakly, up 
to a multiplicative constant, to Y, as N — > oo. 



15 



16 



Step 2. We deal with the case that < A^> < i. Corresponding to (|4.18|) . we assume that 
K < in the following. For every i E {1, • ■ • , a 1 }, we get from Feller [T5][p.281] and Taqqu [2"9"] 
that 

oo 

^fc 2A *- 2 ~-(2A i -l)- 1 j 2A< - 1 ) as j^oo, (4.26) 
fc=j 

since < A D < \. Then by ((4726]) . 

OO 

£fe 2 ^~(7oj 2D - z , as j^oo, (4.27) 
fc =J 

where Co G J 4ut(M ci ) is a diagonal matrix with diagonal entries 

-(2Ai-l)-\--- ,-(2A d -l)- 1 . 
On the other hand, it follows from (|4.19[) that 

j OO 

r(0) + 2^>(z) = -2 £ r(t). (4.28) 

i — 1 i=j-\-X 

By dHHD, 6T27J) and B^51) . 

j 

r(0) + 2 r(i) - ^KCafj 20 - 1 . (4.29) 



i=i 



By (HTJU), P725]) . (|4~2l)| and ([OS]) , we get that as TV -> oo, 

AT 

^at - -2XC , f^j 2 - D - 7 - -2 J fCC* C*ifAr 213 . (4.30) 

3=1 

Since both Co and Ci are diagonal matrices, we choose C G Atii(R d ) such that (C) 2 = CqC\. 
Therefore 

En ~ —2KCN D TN D * C* , (4.31) 

as AT — > oo. Using the same method as the proof of Theorem 12.21 we get that {<3at(£)} converges 
weakly, up to a multiplicative constant, to Y, as N — > oo. The proof has been completed. □ 
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